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Abstract. This paper introduces a simple Markov process inspired by 
the problem of quasicrystal growth. It acts over two-letter words by ran- 
domly performing flips, a local transformation which exchanges two con- 
secutive different letters. More precisely, only the flips which do not in- 
crease the number of pairs of consecutive identical letters are allowed. 
Fixed-points of such a process thus perfectly alternate different letters. 
We show that the expected number of flips to converge towards a fixed- 
point is bounded by 0(n 3 ) in the worst-case and by 0(n 5/2 In n) in the 
average-case, where n denotes the length of the initial word. 



Introduction 

Tilings are often used as a toy model for quasicrystals, with minimal energy 
tilings being characterized by local properties called matching rules. In this con- 
text, a challenging problem is to provide a theory for quasicrystals growth. One 
of the proposed theories relies on a relaxation process ([4] p. 356): a tiling with 
many mismatches is progressively corrected by local transformations called flips. 
Ideally, the tiling eventually satisfies all the matching rules and thus shows a qua- 
sicrystalline structure. It is compatible with experiments, where quasicrystals are 
obtained from a hot melt by a slow cooling during which flips really occur. It 
is however unclear whether only flips can explain successful coolings or if other 
mechanisms should also be taken into account. This question is deeply related 
with the convergence rate of such a flip-correcting process. 

A relaxation process which aims to be physically realist is described in [1]. 
It considers general cut and project tilings of any dimension and codimension, 
and performs each flip which modifies by AE the energy of the tiling with a 
probability depending not only on AE but also on a temperature parameter T 
such that the stationary distribution (at fixed T) is the Boltzmann distribution. 
In this paper, we focus on a very rough version of this general process. First, 
we consider only tilings of dimension and codimension one, which correspond to 
two-letter words. Second, the flips with a corresponding AE being above a fixed 
threshold are equiprobably performed, while the other flips are simply forbidden. 
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The paper is organized as follows. In Sec. 1, we clearly and formally state 
the definition of the stochastic process we consider, as well as the main question 
we are interested in, that is its convergence time. In Sec. 2, we state and prove 
the main result of this paper, which shows that the expected number of flips 
performed by the process to converge is cubic in the length of the initial word 
(Theorem 1). The proof mainly relies on a well-chosen function, called variant, 
whose expected value strictly decreases. This result is then extended in Sec. 3 to 
the case where the initial words are randomly chosen according to some partic- 
ular distribution which aims to be physically realist (Theorem 2). We conclude 
the paper by a short section discussing perspectives. 

1 Settings 

A configuration of length n is a word w = w\ . . . w n over {1, 2} with = |w|2, 
i.e., with as many occurrences of the letter 1 as of the letter 2. Such an object 
is also sometimes called a grand Dyck path, or a bridge. We denote by W n the 
set of configurations of length n. 

It is convenient to represent a configuration w as the broken line of the Eu- 
clidean plane linking the points {k, \wi ■ ■ ■ Wk\i — \wi ■ ■ ■ Wk\2)k=o...,\w\ ( see Fig- !)• 




Fig. 1. The configuration 11211121222112212222211222112111211212. 



A flip at position 1 < i < n on a configuration w e W„ is the local transfor- 
mation which exchanges Wi and Wi + i, provided that these letters are different. 
The height of such a flip is the integer \wi ■ ■ ■ iUj_i|i — |u>i • • • un-\\2- Geome- 
trically, performing a flip on a configuration corresponds to move upwards or 
downwards by 2 a point of the broken line which represents the configuration. 

Flip thus acts over configurations, and it is not hard to check that any two 
configurations with the same length are connected by a sequence of flips. 

A configuration w is said to have a mismatch at position i if Wi = lUj+i. For 
example, the configuration on Fig. 1 has 18 mismatches, while configurations on 
Fig. 2 both have 4 mismatches. The total number of mismatches of w is denoted 



Fig. 2. A flip of height 2 performed at position 5 (leftwards). On the representation, 
the height and the position of a flip are the abscissa and the ordinate of the center 
between the two positions of the moved vertex (that is, the middle of the arrow). 

by E(w). One easily checks that it ranges from for configurations (12)™ and 
(21) n to 2(n - 1) for configurations l n 2" and 2™1". Physically, E(w) can be 
though as the energy of the configuration w, with the configurations (12)" and 
(21) ra thus being ground states. 

We now define a Markov chain on W„ that we call cooling process. It starts 
from wq £ W n and produces a sequence (u>t)t=i,2,... defined by: 

= (w t i£ E(w t ) = 0, 
* +1 \ w' otherwise, 

where w' is obtained by performing on wt a flip uniformly chosen among the 
flips which do not increase the number of mismatches 3 . Note that a flip modifies 
the number of mismatches by at most 2. Flips which decrease the number of 
mismatches are said to be irreversible (they cannot be performed back) while 
the ones which do not modify the number of mismatches are said to be reversible. 




Fig. 3. Flips can increase the number of mismatches (red arrows), let it unchanged (blue 
arrows) or decrease it (green arrows). The two last are respectively called reversible 
and irreversible flips: they are the only ones performed during the cooling process. 

Since the number E(wt) of mismatches is non- increasing, it is natural to ask 
whether it reaches zero or not, and at which rate. We therefore introduce the 
convergence time of this process, which is the random variable T counting the 



3 In other words, the AE threshold discussed into the introduction is equal to zero. 
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number of steps required to transform the initial configuration w into one of 
the two stable configurations without any mismatch, that is 

T := min{t > | E{w t ) = 0}. 

Then, the worst expected convergence time is defined by 

f(n) := max E(T\w = w). 

We are not only interested in the worst case, but also in the average case. Given 
a distribution ^ over W„, we define the \i-averaged expected convergence time by 

T M ■= J2 a*H e ( t K = w). 

tuew„ 



2 Bounding the worst expected convergence time 

In this section, we prove the main result of this paper: 
Theorem 1. The worst expected convergence time is cubic: 

f(n) = 0(n 3 ). 

We will rely on the following proposition (proven, e.g., in [2]) which bounds 
from above the expected convergence time of so-called variants: 

Proposition 1. Consider a Markov chain (x t )t>o over a space ft and a positive 
real map <f> : ft — >• [a, b], called variant. Assume that there is e > such that 
whenever 4>{x t ) > a, E[0(x t+ i) — (j)(x t )\xt\ < —e. Then, the expected value of the 
random variable T := min{t|</>(x 4 ) = a} satisfies 

E(T) < 

e 

Variants are thus sort of "strict" supermartingales, and we need to define a 
suitable one in order to bound the expected convergence time towards mismatch- 
free configurations. Since the number E of mismatches is non-increasing by def- 
inition of the cooling process, we are naturally tempted to use it as a variant. 
However, there are configurations having only reversible flips, whose expected 
variation of mismatches is thus equal to zero (as, for example, w = 1211212212). 
In order to refine, we introduce the notion of Dyck factors (see Fig. 4): 

Definition 1. Consider a configuration w = p-v s. The factor v = V\ ■ ■ ■ Vk is 
said to be a positive Dyck factor of length k and height \p\i — \p\2 if: 

\v\i = \v\ 2 , Vi e {l,...,k}, \v!---Vi\i > \v 1 ---v i \ 2 , - \p\ 2 > 0. 

Exchanging letters 1 and 2 in the above conditions defines a negative Dyck factor. 
A Dyck factor is maximal if no Dyck factor with the same height contains it. 
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Fig. 4. A positive maximal Dyck factor of height 2 and length 6 (in blue), and a 
negative maximal Dyck factor of height —1 and length 12 (in red). 

Let us stress that Dyck factors are either positive or negative, i.e., their 
representations do not cross the y = axis. In particular, any configuration can 
be decomposed in its maximal Dyck factors of height zero, with signs alternating. 

We use Dyck factors to define our variant, which depends on a parameter a: 
Definition 2. For a £ (0, 1), let <j) a be defined on a configuration w by: 

<t> a {w) — J2 ( 1 + Mir, 

v£DF{w) 

where DF(w) denotes the set of maximal Dyck factors ofw. 

Contrary to the number E of mismatches, (j> a can be increased by some flips. 
However, the concavity of x — > x a , which gives to small Dyck factors a weight 
proportionally bigger than the weight given to long Dyck factors, will ensure 
that the average variation over all the flips of a configuration is always negative. 
More precisely: 

Lemma 1. Let w € W^n- If E(w t ) > 0, then the variant 4> a satisfies 

E(0 Q ( W4+1 ) - Mm)\m) < - a{1 2 a) n a ~ 2 . 

Proof. Consider a configuration w = ui\- ■ ■W2n, with E(w) > 0. This ensures 
that at least one flip (reversible or irreversible) can be performed on w. 

Let us first assume that only reversible flips can be performed on w, that is, 
neither 1122 nor 2211 are factors of w. This is equivalent to say that maximal 
Dyck factors contain at least four letters. We will group these reversible flips by 
pairs and prove that the average variation of 4> a over each pair is bounded from 
above by ^n"" 2 . 

Each flip of positive height which transforms 21 into 12 (it increases <f> a ) is 
between two positive flips which transform 12 into 21 (they decrease <j> a ). Re- 
versibility ensures that one of these two flips has the same high as the central 
flip, while the other one is higher; we group the central flip and its higher neigh- 
bor. We proceed symetrically for flips of negative height. Fig. 5 illustrates this. 



Fig. 5. Each flip which increases cf> a (red arrows) can be paired (light green boxes) 
with a flip which decreases <j> a (blue arrows), such that the average variation of <j> a over 
the two flips is negative. Isolated flips can only decrease <f> a . 



Consider two paired reversible flips, say of positive height (the negative height 
case is symmetric). Performing the lowest flip increases by 1 the number p of l's 
in the maximal positive Dyck factor of w starting at the first letter between the 
pair of flips, while performing the highest flip decreases by 1 the number q of l's 
in some other maximal positive Dyck factor of w; one moreover has p > q since 
the latter Dyck factor is higher than the former one (see Fig. 6). 
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Fig. 6. Consider two paired flips (central configuration): the maximal Dyck factor 
extended by the flip increasing <f> a (red arrow, performed leftwards) is greater than the 
one shortened by the flip decreasing (j> a (blue arrow, performed rightwards), i.e., p > q. 



Thus, the average variation of <f> a is 

A-.= 1 -{(jp+2r -(p+ir +q a -(q+m- 

Since x — > x a is concave for a € (0, 1), p > q yields 



a<a :=-((p+2r-(p+ir+p a -(p+ir). 

For k > 1, let Uk — a ( a - 1 )—( a ~ fc + 1 ) _ Note that U2k < 0. One computes: 



2A> / 1 



k>l 



P + l, 

Uk 

(p+iy 



p + i 

= i + y Uk +i + V(-i) fc Uk 



E2 "2fc 
(p + 1 

k>l ^ 



< 2 



k>l 

U2 a(a — 1) 



{ (P+l) 2fc " (P+1) 2 (P + 1) 



2 ' 
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With p < n, this yields the claimed bound for the pair: 

A < A < < (P + X )° x "("- l) = .^-^(p + < _«( 1 -«) w q-a. 

2 (p + 1) 2 2 ^ ; 2 

Last, note that there are some flips which have not been paired, namely the ones 
which are immediatly after or before a crossing of the y — axis. They however 
do not cause trouble, because each of them decreases a maximal Dyck factor of 
length at most n, hence decreases <f> a by at least 

(n - l) a -n a < -an - 1 < - " (1 ~ a) n a ~ 2 . 
v ) 2 

Let us now consider the general case, i.e., when there are also irreversible flips. 
We will show that an irreversible flip increases <p a lesser than (resp. decreases <p a 
more than) two reversible flips do. Hence, the average variation of <j) a can only 
be smaller than in the previous case (which is already negative). Intuitively, it 
is as if each irreversible flip is splitted into two reversible flips before performing 
the above pairing process. 

Let us be more precise. Consider, first, the case of a positive flip which transforms 
21 into 12. It replaces two maximal Dyck factors, say with respectively p and q 
letters 1, by one with p + q + 1 letters 1 (see Fig. 7). Thus, on the one hand: 

A = ( P+q +2r-(p+ir-( q + i) a . 

On the other hand, the variation that would result of two reversible flips is: 

A' = (p + 2) a - [p + l) a + (q + 2) Q - (g + 1)". 
One thus has A < A', since for < a < 1: 

(p + 2)" + (q + 2) a > ((p + 2) + (q + 2)) a = (p + q + 4) Q > (p + q + 2) a . 




Fig. 7. An irreversible flip which replaces two maximal Dyck factors with respectively 
p and q letters 1 by a single maximal Dyck factor with p + q + 1 letters 1. 



We then consider the case of a positive irreversible flip which transforms 12 
into 21. This just removes a maximal Dyck factor with one letter 1 (see Fig. 8). 
Thus, on the one hand: 

A=-2 a . 



Fig. 8. An irreversible flip which removes a maximal Dyck factor with one letter f . 



On the other hand, each of the two neighbor flips of this flip (one on the left, 
the other on the right, both transforming 21 into 12) either decreases the variant 
if it is negative, or increases it by: 

A' =(p+ 2) a - (p + 1) Q < 3" - 2". 

The average variation of 4> a is thus bounded from above by: 

2A' + A 2(3 Q - 2 Q ) - 2" 
3 ~ 3 ' 

A computation shows that, for n > 2, this quantity is lesser than — n°~ 2 
for any a <E (0, 1), with equality for a = 1. The case of a negative irreversible 
flip is symmetric. This concludes the proof. □ 

This lemma thus provides a suitable e (depending on a and n) for Prop. 1. Let 
us now bound <p a and show that it is minimal for mismatch-free configurations: 

Lemma 2. For any configuration w G Wam one has 

(l + n) a <<j> a (w) < (2n) a +\ 

with the lower bound beeing reached if and only if E(w) = 0. 

Proof. Let us first focus on the upper bound. We prove by induction on \ v\ that, 
for any positive configuration v, 

(t> a (v) < (i + Mir +1 - 

The inequality is satisfied for n = 2 since one has <f> a (12) = 2 a . Let us now 
consider a configuration v of length n, and assume that the result holds for any 
shorter positive configuration. Let ui,...,u p be the maximal (positive) Dyck 
factors of v of height 1. One thus has 

v v 

<j> a {v) = (i + MO" + ^) ^ i 1 + Mi)" + E( x + Ni) Q+1 - 

We then use the convexity of x — > x a+1 and compute 



Mv) < (i + H0 a + 2^ 1 + l^li =0- + + Mi < (i + \vUY 
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The claimed result follows by induction. It also holds for negative configurations. 
Let us now consider the case of a general configuration w. Let v\,...,v q be 
its maximal Dyck factors of height 0. Since they are alternatively positive and 
negative configurations, one relies one the previous result to compute 

q q / q \ "+ 1 

<t> a (w) = J2 MVi) < ^(1 + |«,|l) Q+1 < E 1 + = (« + l W ll) a+1 - 

i=l i=l \i=l / 

Since q < \1v\2, this yields 4> a {w) < |w| Q+1 = (2n) a+1 , hence the upper bound. 

Let us now focus on the lower bound. Let robea configuration and vi,...,v g 
be its maximal Dyck factors of height 0. By relying on the concavity of x — > x a , 
one computes 

q q / q \<* 

0aM = E>a(«i) = E;(i + Mi) a > E 1 + Mi =(9 + Hi) a - 

i=l i=l \i=l / 

Since \w\\ — n and q>l, this proves the claimed lower bound. Moreover, if this 
bound is reached, then q = 1, that is, w is either positive or negative. In both 
cases, w has a maximal Dyck factor with n letters 1, which has weight (1 + n) a . 
This ensures that there is no other maximal Dyck factors. Hence, w is one of the 
configuration (12)™ or (21)™, which are exactly the mismatch-free ones and are 
both mapped onto (1 + n) a by <fi a . This proves the claimed result. □ 

This bound is optimal up to a factor 2 Q+1 since one computes 4> a (l n 2 n ) = n a+1 . 
Th. 1 then follows from Prop. 1 with e = ^^-n a - 2 and f = 4> a {w ) < n a+1 : 

E(f (n)) < n^ 1 x = 

a(l — a) a(l — a) 

Here, the best value for a is \ . Other values are useful only in the next section. 



3 Bounding the weighted expected convergence time 

Given a distribution [1 over W„, one deduces from Prop. 1 a bound on the 
^-weighted expected convergence time, with e = "^'^ n"' 2 (Lem. 1): 

tm < I E mM<^M- 

i»ew„ 

Thus, to bound T^(n), we just need to compute the yu-weighted value of cp a . It 
is however generally easier to compute the ^(-weighted value of the volume. The 
volume V(w) of a configuration w is the area between its representation and the 
horizontal axis (for example, Fig. 2, both configurations have volume 15). One 
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can then deduce bound on the /z-weighted value of <j> a from the simple following 
inequality, which holds for any a G (0, 1) and for any configuration w: 

<t> a (w) < V(w). 

For the uniform distribution v over W n , the u-weighted value of the volume 
is known to be equivalent to ^Ln-^/n (see, e.g., [3], p. 533). This allows to prove: 

Theorem 2. The uniformly-weighted expected convergence time T v satisfies 

T v (n)=0(n 5/2 Inn). 

Proof. For any a € (0, 1), one has 



T v (n)<^^x V v(w)V(w) = A y />-«. 

One checks that this upper bound, seen as a function of a, is minimal in 1 — 
For such an a, one has n 7 / 2 ~ a = n 1 / 1 "™^ 5 / 2 — en 5 / 2 . This yields 



T v ( n ) < en 5 / 2 - — ^^- n 5 / 2 Inn = 0(n 5 / 2 Inn). 

V In n ' In n V Inn' 

This shows the claimed result. □ 

However, recall from the introduction that the context of this problem is 
the cooling of a quasicrystal from melt at very high temperature. In particu- 
lar, the natural distributions over initial configurations should be the one of the 
melt. In the melt, the very high temperature imposes a very low threshold for 
AE, so that all the flips are equiproblable, with no flip being any more forbidden. 

Hence, we define the natural distribution v over W n as the stationary dis- 
tribution of the process which, at each step, performs uniformly at random one 
of the possible flips, without any restriction (contrarily to the cooling process 
described in Sec. 1). It is a classical result of random walks on graph that this 
stationary distribution gives to each configuration a weight proportional to the 
number of flips which can be performed onto. In particular, the configurations 
of maximal weight are (12)" and (21)™, i.e., the stable configurations of the 
cooling process, while the configurations of minimal weight are l"2 n and 2™1", 
i.e., the worst configurations of the cooling process. One can thus hope that the 
naturally-weighted expected convergence time T v is lower than the uniformly- 
weighted one. But we only get a similar bound, because the naturally-weighted 
value of the volume has the same growth order as its uniformly- weighted value: 

Proposition 2. The naturally-weighted value of the volume satisfies 



E 



v(w)V(w) ~ — - — n\fn. 
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Proof. We first introduce a combinatorial enumerative series / describing only 
the positive configurations. Let a n ,p, g denotes the number of positive configura- 
tions of size 2n, volume q and onto which p flips can be performed. We set: 

f(z,u,v)= ^2 a n ^ q z n u p v q . 

(n,p, 9 )eN 3 

Let us find a functional equation for /. Each positive configuration w can be 
written w — lu2v, where both u and v are positive configurations, possibly 
empty. The volume satisfies V(w) = V(u) + \u\ + 1 + V(v). The number F of 
flips satisfies either F(w) = F(u) if \v\ = 0, or F(uu) = F(u) + 1 + F(v) otherwise. 
This leads, using symbolic methods described in [3], to the functional equation 

f(z, u, v) = (zuv + zvf(zv, u, v)){l + uf(z, u, v)). (1) 

We now introduce a similar combinatorial enumerative series g describing all the 
configurations. Each configuration w can be written w = v\ ■ ■ ■ v p , where the Vi's 
are alternating positive or negative configurations, hence described by /. Once 
again, symbolic methods described in [3] thus lead to the functional equation 

2z(vf(zv,u,v) + uv) 

g(z,u,v) = - -, 5 —. r r. (2) 

1 — {zuv + zu z v + zvj(zv, U, V) + zuv) 



The naturally- weighted value of the volume can then be derived from g as follows: 



f \-\rl \ L " ' dvdu9 V" 2 ' U ' V ) \u=l,v=l 

^ v(w)V(w) = ° ° . 



where [z n ]h denotes the coefficient of z n in the series h. Let us compute this. 

We first consider the denominator. The differentiation of (2) with respect to 
u yields for || (z, u, v) a huge formula where appear f(zv, u, v) and ^f(zv, u, v). 
On the one hand, f(z, 1, 1) is simply the usual series describing Dyck paths: 



/(*,M)= Yz • 

On the other hand, the differentiation of (1) with respect to u yields 

df ( df \ 

— (2;, u,v) = izv + vz-^-(zv, u, v)\(l + uf(z, u, v)) 

( df 

+ (zuv + zvf(zv,u, v)) I f(z,u, v) + u—(z,u,v) 

Since f(z, 1, 1) is known, this allows to compute 



df df 1 - y/T=te - 2Z 
— (Z, U,V) „=!,„=! = -^-(zV,U,v)\ u=hv=1 = . 

ou du 2zy/l - Az 
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Working this out in the equation for §§(z, u, v), one gets 



For the numerator, a similar but even more tedious computation yields 
AT , . d 2 g . .. 16z 3 + 4z 2 + 2z 

N{z) := 9^; (z ' w ' v)l - i -= 1 = (i - 4z)3 • 

Both D(z) and N(z) are Z\-analytic at z = \. In this point, one computes 
D{z) ~ — \ . , /2 and iV(z) 



2(!_4 z )3/2 — w ( 1 _ 4z )3- 

Asymptotics of the coefficient of z n in the series expansions of D and iV can 
then be deduced from the following result (see e.g. [3], Chap. VI), which holds 
for a £ {-1, -2, . . .} and p <= C: 

l^'-f-fS'"'"-' 

Indeed, with p equals to \ and a respectively equals to | and 3, this yields: 

[zn]D{z) ~ 27W2) 4nV " and [z " ]iV(z) ~ rk 4 "" 2 ' 

One thus finally computes 

, .... , [z n ]N(z) 2r(3/2) 

tuGW 2 „ w 

The claimed result follows by replacing 2n by n (the sum is taken over W„) □ 



4 Perspectives 

Th. 1 and 2 provide wpper bounds on the expected convergence time of the cool- 
ing process, respectively in the worst case and in the average case (uniform and 
natural distributions). Experiments suggest, especially in the worst case, that 
these bounds are tight (see Fig. 9 and 10). This remains to be rigorously proven. 

Moreover, this paper focuses only on the expectation of the convergence time, 
as the physically most significant moment. It would be worth to consider also 
higher moments, in order to obtain the limit distribution of the convergence time. 

Last but not least, twodetter words correspond to tilings of dimension one 
and codimension one. But the physically interesting cases have dimension two or 
three and codimension at least two. For example, the celebrated Penrose tiling 
has dimension two and codimension three. The analysis of these cases seems 
harder but rewarding: experiments show much faster convergence rates (see [1]). 
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Fig. 9. The blue points represent the average over 10 simulations of the worst conver- 
gence time (on the x-axis, the length n of configurations). The pink curve 0.17n 3 fits 
remarkably well with experimental values. The standard deviation of the worst conver- 
gence time is experimentally very small, so that only 10 simulations already provide a 
good approximation of its expectation. 
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Fig. 10. The blue points represent the average over 1000 simulations of the uniformly- 
weighted convergence time (on the x-axis, the length n of configurations). The pink 
curve 0.24n 5 ^ 2 Inn fits rather well with experimental values. The standard deviation 
of the uniformly-weighted convergence time is experimentally much greater than for 
the worst convergence time (indeed, each simulation starts from a configuration cho- 
sen uniformly at random), so that the number of simulations needed to get a good 
approximation is much greater too. 
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